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Abstract 

We use the coupled cluster method (CCM) to study the zero-temperature phase diagram of a 
two-dimensional frustrated spin-half antiferromagnet, the so-called Union Jack model. It is defined 
on a square lattice such that all nearest-neighbor pairs are connected by bonds with a strength 
Ji > 0, but only half the next-nearest-neighbor pairs are connected by bonds with a strength 
J2 = nJi > 0. The bonds are arranged such that on the 2x2 unit cell they form the pattern of 
the Union Jack flag. Alternating sites on the square lattice are thus 4-connected and 8-connected. 
We find strong evidence for a first phase transition between a Neel antiferromagnetic phase and a 
canted ferrimagnetic phase at a critical coupling = 0.66 it 0.02. The transition is an interesting 
one, at which the energy and its first derivative seem continuous, thus providing a typical scenario of 
a second-order transition (just as in the classical case for the model), although a weakly first-order 
transition cannot be excluded. By contrast, the average on-site magnetization approaches a nonzero 
value = 0.195 it 0.005 on both sides of the transition, which is more typical of a first-order 
transition. The slope, dM/dn, of the order parameter curve as a function of the coupling strength 
K, also appears to be continuous, or very nearly so, at the critical point thereby providing 
further evidence of the subtle nature of the transition between the Neel and canted phases. Our 
CCM calculations provide strong evidence that the canted ferrimagnetic phase becomes unstable 
at large values of k, and hence we have also used the CCM with a model collinear semi-stripe- 
ordered ferrimagnetic state in which alternating rows (and columns) are ferromagnetically and 
antiferromagnetically ordered, and in which the spins connected by J2-honds are antiparallel to 
one another. We find tentative evidence, based on the relative energies of the two states, for a 
second zero-temperature phase transition between the canted and semi-stripe-ordered ferrimagnetic 
states at a large value of the coupling parameter around ~ 125 ±5. This prediction, however, is 
based on an extrapolation of the CCM results for the canted state into regimes where the solutions 
have already become unstable and the CCM equations based on the canted state at any level of 
approximation beyond the lowest have no solutions. Our prediction for is hence less reliable 
than that for Kcj- Nevertheless, if this second transition at does exist, our results clearly 
indicate it to be of first-order type. 

PACS numbers: 75.10.Jni, 75.30.Gw, 75.40.-s, 75.50.Ee 
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I. INTRODUCTION 



Quantum magnetism at zero temperature for lattices in two spatial dimensionsi"- is an 
important and fascinating subject because such systems display a wide variety of behavior, 
including semi-classical Neel ordering, two-dimensional quantum "spirals", valence-bond 
crystals/solids, and spin liquids. The behavior of these systems is driven by the nature of 
the underlying crystallographic lattice, the number and range of bonds on this lattice, and 
the spin quantum numbers of the atoms localised to the sites on the lattice. There are 
very few exact results for quantum spin systems on two-dimensional (2D) lattices, and so 
the application of approximate methods is crucial to their understanding. The theoretical 
investigation of these models has been strongly mirrored by the discovery and experimental 
investigation of new quasi-2D magnetic materials. It seems clear that we can only form a 
complete picture of such 2D quantum spin-lattice systems by considering a wide range of 
possible scenarios that are often inspired (or followed shortly afterwards) by experimental 
studies. 

A prototypical case is presented by the spin-half square-lattice Heisenberg antiferromag- 
net (HAF) model. This model has been studied extensively via a range of approximate 
techniques.-^- Its basic properties have been well-established, where, for example, approxi- 
mate results for the order parameter indicate that about 61% of the classical Neel ordering 
persists in the quantum limit at zero temperature. A review of the properties of the spin- 
half square-lattice HAF is given in Ref. j7|. The most accurate results for this model are 
provided by quantum Monte Carlo (QMC) simulations.- Indeed, QMC techniques generally 
provide the benchmark for quantum magnets in two spatial dimensions. However, its use is 
severely limited by the "sign problem," which is often a consequence of quantum frustration 
in the context of lattice spin systems. 

A common theme has also begun to emerge recently when frustrating next-nearest- 
neighbor (NNN) bonds with strength J2 > are added to the basic spin-half square-lattice 
HAF with nearest-neighbor (NN) bonds with strength Ji > 0. The frustrating J2 bonds may 
be added on some or all of the square plaquettes of the lattice and/or across both or only 
one of the diagonals of each plaquette. Perhaps the prototypical such model is the so-called 
spin-half J1-J2 model in which all possible NNN bonds are included. Recent interest in this 
model has been reinvigorated by the discovery of various layered magnetic materials, such 
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as Li2VOSi04, Li2VOGe04, VOM0O4, and BaCdVO(P04)2. Several approximate methods 
have been used to simulate the properties of this system including the coupled cluster method 
(CCM),-~— series expansion (SE) techniques,—"— exact diagonalization (ED) methods,—"— 
and hierarchical mean-field (MF) calculations.— These approximate techniques have estab- 
lished conclusively that there are two phases exhibiting magnetic long-range order (LRO) 
at small and at large values of k = J2/J1 respectively. For k < ~ 0.4 the ground-state 
(gs) phase exhibits (NN) Neel magnetic LRO, whereas for k > ~ 0.6 it exhibits coUinear 
striped LRO in which alternating rows (or columns) of the square lattice have opposite spins, 
with the spins on each row (or columns) aligned, so that the Neel order is between NNN 
pairs. The intermediate region consists of a quantum paramagnetic state without magnetic 
LRO. 

Several other models in this general class of spin-half models with both NN and NNN 
interactions have prompted recent interest. They all involve the removal of some of the NNN 
J2-bonds from the fundamental J1-J2 model. One such example is the Shastry-Sutherland 
model,—"— realized experimentally by the magnetic material SrCu(B03)2, which involves 
the removal of three-quarters of the J2-bonds. Whereas the J1-J2 model on the 2D square 
lattice has each of the sites connected by 8 bonds (4 NN Ji-bonds and 4 NNN J2-bonds) 
to other sites, the Shastry-Sutherland model has each of the sites connected by 5 bonds 
(4 NN Ji-bonds and 1 NNN J2-bond). A second example is the HAF on the anisotropic 
triangular lattice model (also known as the interpolating square-triangle model),— realized 
experimentally by the magnetic material CS2CUCI4, which involves the removal of half the 
J2-bonds from the original J1-J2 model. In this model each of the sites on the 2D square 
lattice is connected by 6 bonds (4 NN Ji-bonds and 2 NNN J2-bonds) to other sites, such that 
the remaining J2-bonds connect equivalent NNN sites in each square plaquette. Although 
all of the models mentioned above show antiferromagnetic Neel ordering for small J2, their 
phase diagrams for larger J2 display a wide variety of behavior, including, two-dimensional 
quantum "spirals", valence-bond crystals/solids, and spin liquids. Thus, in the absence 
of any definitive theoretical argument, the best way to understand this class of NN/NNN 
models on the square lattice is to treat each one on a case-by-case basis. 

In this article we study another frustrated spin-half model that has both NN (Ji) and 
NNN ( J2) bonds on the square lattice, where these bonds form a pattern that resembles the 
"Union Jack" fiag. Just as for the anisotropic triangular HAF described above, the Union 
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Jack model on the 2D square lattice also has only one frustrating NNN bond per square 
plaquette, but these J2-bonds are now arranged such that half the sites are connected by 8 
bonds (4 NN Ji-bonds and 4 NNN J2-bonds) to other sites, while the other half are connected 
only by 4 Ji-bonds to their NN sites, as described more fully below in Sec. [TTl This model has 
previously been studied using SWT— and SE techniques.— As in the case of the spin-half 
interpolating square-triangle model, it was shown^^"— that NN Neel order for the Union Jack 
model persists until a critical value of the frustrating NNN ( J2) bonds. However, in contrast 
to the interpolating square-triangle model, there exists a ferrimagnetic ground state in which 
spins on the eight-connected sites cant at a nonzero angle with respect to their directions 
in the corresponding Neel state. This model thus exhibits an overall magnetic moment in 
this regime, which is quite unusual for spin-half 2D materials with only Heisenberg bonds 
and which therefore preserve (spin) rotational symmetries in the Hamiltonian. This model 
also presents us with a difficult computational task in order to simulate its properties. Here 
we wish to study this model using the CCM, which has consistently been shown to yield 
insight into a wide range of problems in quantum magnetism, and which we now hope will 
hence shed yet more light on the whole class of NN/NNN models as mentioned before. The 
dual associated features of a model with two sorts of sites with differing connectivities, and 
its consequent ferrimagnetic phase, are just those that are attracting the interest of the 
community now. 

II. THE MODEL 

In this paper we now apply the CCM to the spin-half Union Jack model that has been 
studied recently by other means.—"— Its Hamiltonian is written as 



where the operators Sj = (sf,s^,s^) are the quantum spin operators on lattice site i with 
sf = s{s + 1) and s = 1/2. On the square lattice the sum over runs over all distinct 
NN bonds with strength Ji, while the sum over [i, k] runs over only half of the distinct NNN 
diagonal bonds having strength J2 and with only one diagonal bond on each square plaquette 
as arranged in the pattern shown explicitly in Fig. [H The unit cell is thus the 2x2 square 
shown in Fig. [Tal (We note that, by contrast, the J1-J2 niodel discussed above, includes all 
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(a) Canted (b) Scmi-stripcd 

FIG. 1: (Color online) Union Jack model; — Ji; J2. (a) Canted state; (b) 

Semi-striped state. The unit cell is a square of side length 2. 

of the diagonal NNN bonds on the square lattice.) We consider here the case where both 
sorts of bonds are antiferromagnetic, Ji > and J2 = nJi > 0, and are hence acting to 
compete against (or to frustrate) each other. Henceforth we set Ji = 1. We consider the 
model equivalently defined by the Union Jack geometry in which there are two sorts of sites, 
namely the A sites with 8 NN sites and the B sites with 4 NN sites, as shown in Fig. [Tal 

Considered classically rather than quantum-mechanically, (and thus corresponsing to the 
quantum case in the limit where the spin quantum number s — ?■ 00), the Union Jack model 
has only two gs phases as the parameter k is varied over the range (0, 00). A simple 
variational analysis for the classical model reveals that for < k < 1/2 the gs phase is Neel- 
ordered, exactly as for the full J1-J2 model. Thus the Neel ordering induced by the Ji-bonds 
acting alone is preserved as the strength of the competing J2-bonds is increased, until the 
critical value = 0.5 is reached. For k > nf a new phase of lower energy emerges, just as in 
the full J1-J2 niodel. However, whereas for the full J1-J2 model that new phase is a classical 
striped state in which alternate rows (or columns) of spins are arranged antiparallel to one 
another, the new classical gs phase for the Union Jack model is the canted ferrimagnetic 
state shown in Fig. [la] in which the spins on each of the alternating Ai and A2 sites of the A- 
sublattice are canted respectively at angles (7r=i=0) with respect to those on the B-sublattice, 
all of the latter of which point in the same direction. On the A-sublattice each site Ai has 
4 NN sites A2, and vice versa. The angle between the NN spins on the A-sublattice is thus 
2(f). 
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The classical energy of the above canted state is thus 

E = A^s^(/tcos20 - 2 cos 0) , 



(2) 



where Ji = 1 and — )■ oo is the number of sites. Clearly the energy is extremized when 

sin0(l — 2fi;cos0) = . (3) 

When K < nf = 0.5, the lowest energy corresponds to sin0 = and hence to the Neel state. 
By contrast, when k, > K,f = 0.5 the lowest energy solution is the canted state with 



cos- ( i- ) . (4) 



Thus the classical gs energy is given by 



E'^ = I ^ ' " (5) 

The classical phase transition at k = nf = 0.5 is of continuous (second-order) type with the 
gs energy and its first derivative both continuous functions of k, although there are finite 
discontinuities in the second- and higher-order derivatives at k = nf. 

In the classical canted phase the total magnetization per site is m^^ = — (2fi;)~^], 
and the model thus exhibits ferrimagnetism in this phase. Whereas ferrimagnetism more 
commonly occurs when the individual ionic spins have different magnitudes on different 
sublattices, it arises here in a case where when the spins all have the same magnitude and 
all the interactions are antiferromagnetic in nature, but the frustration between them acts to 
produce an overall magnetization. The total magnetization m vanishes linearly as k — 
from the canted phase and then remains zero in the Neel phase for k, < nf. The spontaneous 
breaking of the spin rotation symmetry is also reflected by the vanishing of the energy gap on 
both sides of the transition. Clearly on both sides of the transition the translation symmetry 
of the lattice is also broken. 

One of the aims of the present paper is to give a fully microscopic analysis of the Union 
Jack model for the quantum case where the spins all have spin quantum number s = 1/2. 
We are interested to map out the zero-temperature (T = 0) phase diagram of the model, in- 
cluding the positions and orders of any quantum phase transitions that emerge. In particular 
we investigate the quantum analogs of the classical Neel and canted phases and calculate the 



effect of quantum ffuctuations on the position and nature of the transition between them. 
We also aim to investigate, for particular regions of the control parameter k, whether the 
quantum fluctuations may favor other phases, which have no classical counterparts. One 
such possible candidate is discussed below. 

In the limit of k — t- oo the above classical limit corresponds to a canting angle — 
^TT, such that the spins on the A-sublattice become Neel-ordered, as is expected. The 
spins on the antiferromagnetically-ordered A-sublattice are orientated at 90° to those on the 
ferromagnetically-ordered B-sublattice in this limit. In reality, of course, there is complete 
degeneracy at the classical level in this limit between all states for which the relative ordering 
directions for spins on the A- and B- sublattices are arbitrary. Clearly the exact spin- 
1/2 limit should also comprise decoupled antiferromagnetic and ferromagnetic sublattices. 
However, one might now expect that this degeneracy in the relative spin orientations between 
the two sublattices is lifted by quantum fluctuations by the well-known phenomenon of order 
by disorder.— Just such a phase is known to exist in the full spin-1/2 J1-J2 model for values 
of J2/J1 ^ 0.6, where it is the so-called collinear striped phase in which, on the square 
lattice, spins along (say) the rows in Fig. [1] order ferromagnetically while spins along the 
columns and diagonals order antiferromagnetically. We have also shown how such a striped 
state is stabilized by quantum fluctuations for values of J2/J1 > 1.8 for the spin-1/2 J1-J2 
model defined on an anistropic 2D lattice,— as discussed in Sec. [T] above. 

The existence of the striped state as a stable phase for large values of the frustration 
parameter for both the spin-1/2 J1-J2 and J1-J2 models above is a reflection of the well- 
known fact that quantum fluctuations favor collinear ordering. In both cases the order-by- 
disorder mechanism favors the collinear state from the otherwise infinitely degenerate set of 
available states at the classical level. For the present Union Jack model the corresponding 
collinear state that might perhaps be favored by the order by disorder mechanism is the 
so-called semi-striped state shown in Fig. [lb] where the A-sublattice is now Neel-ordered 
in the same direction as the B-sublattice is ferromagnetically ordered. Alternate rows (or 
columns) are thus ferromagnetically and antiferromagnetically ordered in the same direction. 
We investigate the possibility below that such a semi-stripe-ordered phase may be stabilized 
by quantum fluctuations at larger values of k. 
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III. THE COUPLED CLUSTER METHOD 



The CCM (see, e.g., Refs. 30|-|32| and references cited therein) that we employ here is one 
of the most powerful and most versatile modern techniques available to us in quantum many- 
body theory. It has been applied very successfully to various quantum magnets (see Refs. |8|- 
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23 



25 



32 



and references cited therein). The method is particularly appropriate for 
studying frustrated systems, for which some of the main alternative methods either cannot 
be applied or are sometimes only of limited usefulness, as explained below. For example, 
QMC techniques are particularly plagued by the sign problem for such systems, and the ED 
method is restricted in practice by available computational power, particularly for s > 1/2, 
to such small lattices that it is often insensitive to the details of any subtle phase order 
present. 

The method of applying the CCM to quantum magnets has been described in detail 
elsewhere (see, e.g., Refs. 30|-|35| and references cited therein). It relies on building multispin 
correlations on top of a chosen gs model state |$) in a systematic hierarchy of LSUBn 
approximations (described below) for the correlation operators S and S that parametrize 
the exact gs ket and bra wave functions of the system respectively as |\E') = e'^l'l') and 
= {^\Se~^. In the present case we use three different choices for the model state 1$), 
namely either of the classical Neel and canted states, as well as the semi-striped state. Note 
that for the canted phase we perform calculations for arbitrary canting angle [as shown in 
Fig. [la], and then minimize the corresponding LSUBn approximation for the energy with 
respect to 0, -E'LsuBn(0) min ^ = 0LSUBn- Generally (for n > 2) the minimization 
must be carried out computationally in an iterative procedure, and for the highest values 
of n that we use here the use of supercomputing resources was essential. Results for the 
canting angle 0LSUBn will be given later. We choose local spin coordinates on each site in 
each case so that all spins in |$), whatever the choice, point in the negative z-direction (i.e., 
downwards) by definition in these local coordinates. 

Then, in the LSUBn approximation all possible multi-spin-flip correlations over different 
locales on the lattice defined by n or fewer contiguous lattice sites are retained. The numbers 
of such distinct (i.e., under the symmetries of the lattice and the model state) fundamental 
configurations of the current model in various LSUBn approximations are shown in Table [H 
We note that the distinct configurations given in Table |T]are defined with respect to the Union 
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TABLE I: Number of fundamental LSUBn configurations (Nf) for tlie semi-striped and 
canted states of the spin-1/2 Union Jack model, based on the Union Jack geometry defined 

in the text. 



Nf 

Method 

semi-striped canted 



LSUB2 


3 


5 


LSUB3 


5 


42 


LSUB4 


41 


199 


LSUB5 


194 


1259 


LSUB6 


1159 


8047 


LSUB7 


6862 


56442 



Jack geometry described in Sec. [TTl in which the B-sublattice sites of Fig. [Ta] are defined to 
have 4 NN sites and the A-sublattice sites are defined to have the 8 NN sites joined to them 
either by Ji- or J2-bonds. If we chose instead to work in the square-lattice geometry every 
site would have 4 NN sites. The coupled sets of equations for these corresponding numbers 
of coefficients in the operators S and 5* are derived using computer algebra^ and then 
solved^ using parallel computing. We note that such CCM calculations using up to about 
10^ fundamental configurations or so have been previously carried out many times using the 
CCCM code^ and heavy parallelization. A significant extra computational burden arises 
here for the canted state due to the need to optimize the quantum canting angle at each 
LSUBra level of approximation as described above. Furthermore, for many model states the 
quantum number = Xlili original global spin-coordinate frame, may be used 

to restrict the numbers of fundamental multi-spin-fiip configurations to those clusters that 
preserve as a good quantum number. This is true for the Neel state where = and for 
the semi-striped state for which = N/A, where is the number of lattice sites. However, 
for the canted model state that symmetry is absent, which largely explains the significantly 
greater number of fundamental configurations shown in Table [T] for the canted state at a 
given LSUBn order. Hence, the maximum LSUBn level that we can reach here for the 
canted state, even with massive parallelization and the use of super computing resources, is 
LSUB7. For example, to obtain a single data point for a given value of k (i.e., for a given 
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value of J2, with Ji = 1) for the canted phase at the LSUB7 level typically required about 
0.3 hours computing time using 600 processors simultaneously. However, for values of n 
near to termination points at which CCM solutions using that model state disappear (as 
described more fully below), the computing time typically increased significantly. 

At each level of approximation we may then calculate a corresponding estimate of the gs 
expectation value of any physical observable such as the energy E and the magnetic order 
parameter, M = — X]ili(^kf l^); defined in the local, rotated spin axes, and which thus 
represents the average on-site magnetization. Note that M is just the usual sub lattice (or 
staggered) magnetization per site for the case of the Neel state as the CCM model state, for 
example. 

It is important to note that we never need to perform any finite-size scaling, since all 
CCM approximations are automatically performed from the outset in the infinite-lattice 
limit, — 00, where is the number of lattice sites. However, we do need as a last step 
to extrapolate to the n — )■ 00 limit in the LSUBra truncation index n. We use here the 
well-tested^^i^ empirical scaling laws 

E/N = flo + ciin~'^ + a2?T.~^ , (6) 

M = feo + hn-^ + hn-"" . (7) 

IV. RESULTS 

We report here on CCM calculations for the present spin-1/2 Union Jack model Hamil- 
tonian of Eq. ([T]) for given parameters (Ji = 1, J2), based respectively on the Neel, canted 
and semi-striped states as CCM model states. Our computational power is such that we 
can perform LSUBn calculations for each model state with n < 7. We note that, as has 
been well documented in the past,— the LSUBn data for both the gs energy per spin E/N 
and the average on-site magnetization M converge differently for the even-n and the odd-n 
sequences, similar to what is frequently observed in perturbation theory.— Since, as a gen- 
eral rule, it is desirable to have at least (n + 1) data points to fit to any fitting formula that 
contains n unknown parameters, we prefer to have at least 4 results for different values of 
the LSUBn truncation index n to fit to Eqs. ([6]) and ([7]). However, for all of our extrapolated 
results below we perform separate extrapolations using the even and odd LSUBn sequences 
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2(|)/it 2(|)/ji 

(a) (b) 

FIG. 2: (Color online) Ground-state energy per spin of the spin-1/2 Union Jack 
Hamiltonian of Eq. ([1]) with Ji = 1, using the LSUB6 approximation of the CCM with the 
canted model state, versus the canting angle 0, for some illustrative values of J2 in the 
range < J2 < 1.5 for Fig. [2a] and 0.65 < J2 < 0.75 for Fig. [2bl For J2 < 0.68 in this 
approximation the minimum is at </> = (Neel order), whereas for J2 > 0.68 the minimum 
occurs at = 0lsub6 7^ 0, indicating a phase transition at J2 ~ 0.68 in this LSUB6 

approximation. 

with n = {2, 4, 6} and n = {3, 5, 7}. 

A. Neel state versus the canted state 

We report first on results obtained using the Neel and canted model states. While classi- 
cally we have a second-order phase transition from Neel order (for k < nf) to canted order 
(for K > hif), where k = J2/J1, at a value hif = 0.5, using the CCM we find strong indica- 
tions of a shift of this critical point to a higher value ~ 0.66 in the spin-1/2 quantum 
case as we explain in detail below. Thus, for example, curves such as those shown in Fig. 
|2] show that the Neel model state (</> = 0) gives the minimum gs energy for all values of 
K, < where k^^ = /^LSUBn g^jgQ dependent on the level of LSUBn approximation, as we 
see clearly in Fig. |3l 

By contrast, for k, > the minimum in the energy is found to occur at a value 7^ 0. 
If we consider the canting angle itself as an order parameter (i.e., = for Neel order 
and 7^ for canted order) a typical scenario for a first-order phase transition would be the 
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(a)n = {2,4,6} (b)n = {3,5,7} 

FIG. 3: (Color online) The angle 0LSUBn that minimizes the energy -ELSUBn(0) of the 
spin- 1/2 Union Jack Hamiltonian of Eq. ([1]) with Ji = 1, in the LSUBn approximations 
with (a) n = {2, 4, 6} and (b) n = {3, 5, 7}, using the canted model state, versus J2. The 
corresponding classical result (pci from Eq. (jlj) is shown for comparison. We find in the 
LSUBn quantum case with n > 2 a weakly first-order phase transition or second-order 
phase transition (e.g., for LSUB6 at J2 ~ 0.680 and LSUB7 at J2 ~ 0.646). By contrast, in 
the classical case there is a second-order phase transition at J2 = 0.5. 

appearance of a two-minimum structure for the gs energy as a function of 0. If we therefore 
admit such a scenario, in the typical case one would expect various special points in the 
transition region, namely the phase transition point itself where the two minima have 
equal depth, plus one or two instability points and where one or other of the minima 
(at = and </> 7^ respectively) disappears. By contrast, a second-order phase transition 
might manifest itself via a one-minimum structure for the gs energy as a function of 0, in 
which the single minimum moves smoothly and continuously from the value = for all 
values of K, < to nonzero value (p ^ for n > Hci- 

We show in Fig. [2]our results for the LSUB6 approximation based on the canted (or Neel) 
state as the CCM model state. Very similar curves occur for other LSUBn approximations. 
A close inspection of curves such as those shown in Fig. [2] for the LSUB6 case shows that 
what happens for this model at this level of approximation is that for n < 0.68 the only 
minimum in the gs energy is at = (Neel order). As this value is approached from below 
the LSUB6 energy curves become extremely flat near = 0, indicating the disappearance 
at = of the second derivative (PE/dcfP' (and possibly also of one or more of the higher 
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derivatives dT'E /dcj)^ with n > 3), as well as of the first derivative dE/d(j). Then, for all 
values K > 0.68 the LSUB6 curves develop a minimum at a value 7^ which is also the 
global minimum. The state for 7^ is thus the quantum analog of the classical canted 
phase. The fact that the antiferromagnetic Neel order survives into the classically unstable 
regime is another example of the well-known phenomenon that quantum fluctuations tend 
to promote coUinear order in rnagnetic spin-lattice systems, as has been observed in many 
other such cases (see e.g., Ref. 34|j40| ) . Thus, this coUinear Neel-ordered state survives into 
a region where classically it becomes unstable with respect to the non-coUinear canted state. 

A close inspection of the curves shown in Fig. [3] for various LSUBn approximation shows 
that the crossover from one minimum (0 = 0, Neel) solution to the other (0 7^ 0, canted) 
appears to be continuous for the odd-ra sequence, thus indicating a second-order transition 
according to the above scenario. By contrast, for the even-n sequence with n > 2 the 
curves in Fig. |3] become very steep in the crossover region just above kLSUBti ^.j^g 
extremely fiat nature of the gs energy curves as a function of in this region, as shown in 
Fig. 121 it is impossible to rule out a small but finite discontinuity in the curves of Fig. 123 
for the even-n LSUBn sequence at k = k^^^^^". However, if the phase transition is, in fact 
first-order, it is certainly only very weakly so according to this criterion. 

Thus, based on the evidence presented so far of the gs energies of the Neel and canted 
phases, it would appear that the transiton at k = between these two phases is either 
second-order, as in the classical phase, or weakly first-order. Such a situation where the 
quantum fluctuations change the nature of a phase transition qualitatively from a classical 
second-order type to a quantum first-order type has also been seen previously in the compa- 
rable spin- 1/2 HAF models that interpolate continuously between the square and triangular 
lattices,— and between the square and honeycomb lattices,— respectively. In the present 
spin-1/2 Union Jack model, however, the CCM gs energy results appear to favor a second- 
order transition, although the extreme insensitivity of the gs energy to the canting angle 
near the crossover region, especially for the even-n LSUBn sequence with > 2, means that 
we cannot rule out a weakly first-order transition. The evidence to date indicates, however, 
that the quantum phase transition at i^c-^ IS cl subtle one. Furthermore, the present spin-1/2 
Union Jack model appears, on the evidence to date, to behave somewhat differently (viz., in 
some senses "more classically") than its corresponding spin-1/2 interpolating square-triangle 
Heisenberg antiferromagnet counterpart.— Further evidence from Fig. [3] appears to back up 
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TABLE II: The critical value K^suBn which the transition between the Neel phase 
= 0) and the canted phase (0 7^ 0) occurs in the LSUBn approximation using the CCM 
with (Neel or) canted state as model state. 



Method K^suBn 
LSUB2 0.740 
LSUB4 0.696 
LSUB6 0.680 

LSUBoo 0.651 ±0.001 . 

LSUBoo 0.676 ±0.004 . 
LSUB3 0.597 
LSUB5 0.630 
LSUB7 0.645 

LSUBoo " 0.681 ±0.001 . 

LSUBoo 0.653 ±0.004 . 



" Based on k^^^^"^ = ao + ain'^, with n = {2, 4, 6} or n = {3, 5, 7} . 
Based on K^suBn = 5^ + bin''^, with n = {2, 4, 6} orn = {3, 5, 7} . 

this observation. Thus, we see from Fig. |3] that the quantum canting angle approaches 
its asymptotic value 7r/2 as k — )■ 00 slightly faster than does the corresponding classical 
value. By contrast, in the case of the spin-1/2 interpolating square-triangle Heisenberg 
antiferro magnet,— the corresponding pitch angle </> of the spiral phase (that is the analog 
of the canted phase for the present model) approaches its similar asymptotic value 7r/2 as 
K — > 00 very much faster than does the classical value. We also discuss this difference more 
fully below, where we find further evidence that quantum fluctuations modify the classical 
behavior of the Union Jack model rather less than they do for its corresponding spin-1/2 
interpolating square-triangle Heisenberg antiferromagnet counterpart. 

We show in Table ITT] the critical values K^suBn which the transition between the Neel 
and canted phases occurs in the various LSUBn approximations shown in Fig. [31 In the 
past we have found that a simple linear extrapolation, /^LSUBn _ _|_ ain~^, yields a good 
fit to such critical points, as seems to be the case here too. The corresponding "LSUBoo" 
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estimates from the LSUBn data in Table ITTl are k^^ = 0.651±0.001 based on n = {2, 4, 6} and 
Kci = 0.681 ± 0.001 based on n = {3, 5, 7} where the quoted errors are simply the standard 
deviations from the two fits. Similar estimates based on an extrapolation ^^suBn _ ^^_|_^^j^-2 
are also shown in Table [ITl for which the standard deviations are clearly greater. The fact 
that the two estimates based on the even-n and odd-n LSUBn sequences differ slightly from 
one another is a reflection of the extreme insensitivity of the gs energy to the canting angle 
(j) near k^^^^"', and the difference between the two estimates is a rough indication of our real 
error bars on k^. We also present other independent estimates of below. 

We note from Fig. [2] that for certain values of J2 with Ji = 1 (or, equivalently, k) CCM 
solutions at a given LSUBn level of approximation (viz., LSUB6 in Fig. [2]) exist only for 
certain ranges of the canting angle (p. For example, for the pure square-lattice HAF (k = 0) 
the CCM LSUB6 solution based on a canted model state only exists for < </> < 0.16l7r. In 
this case, where the Neel solution is the stable ground state, if we attempt to move too far 
away from Neel coUinearity the CCM equations themselves become "unstable" and simply 
do not have a real solution. Similarly, we see from Fig. [2] that for n = 1.5 the CCM LSUB6 
solution exists only for O.SOStt < < 0.4457r. In this case the stable ground state is a canted 
phase, and now if we attempt either to move too close to Neel collinearity or to increase the 
canting angle too close to its asymptotic value of 7r/2, the real solution terminates. 

Such terminations of CCM solutions are very common and are very well documented.— 
In all such cases a termination point always arises due to the solution of the CCM equations 
becoming complex at this point, beyond which there exist two branches of entirely unphysical 
complex conjugate solutions.— In the region where the solution reflecting the true physical 
solution is real there actually also exists another (unstable) real solution. However, only the 
shown branch of these two solutions reflects the true (stable) physical ground state, whereas 
the other branch does not. The physical branch is usually easily identified in practice as 
the one which becomes exact in some known (e.g., perturbative) limit. This physical branch 
then meets the corresponding unphysical branch at some termination point (with infinite 
slope on Fig. [2]) beyond which no real solutions exist. The LSUBn termination points are 
themselves also reflections of the quantum phase transitions in the real system, and may be 
used to estimate the position of the phase boundary,— although we do not do so for this 
first critical point since we have more accurate criteria discussed above as well as below. 

Before doing so, however, we wish to give some further indication of the accuracy of our 
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TABLE III: Ground-state energy per spin E/N and magnetic order parameter M (i.e., the 
average on-site magnetization) for the spin-1/2 square-lattice HAF. We show CCM results 
obtained for the Union Jack model with Ji = 1 and J2 = using the Neel model state in 
various CCM LSUBn approximations defined on the Union Jack geometry described in 
Sec. ini We compare our extrapolated {n — )■ 00) results using Eqs. (jH]) and ([7]) with the 
odd-77, and even-n LSUBn data sets with other calculations. 



Method 


E/N 


M 


LSUB2 


-0.64833 


0.4207 


LSUB3 


-0.65044 


0.4151 


LSUB4 


-0.66366 


0.3821 


LSUB5 


-0.66398 


0.3795 


LSUB6 


-0.66703 


0.3630 


LSUB7 


-0.66724 


0.3606 


Extrapolations 


LSUBoo 


-0.6698 


0.316 


LSUBoo 


^ -0.6704 


0.304 


QMC " 


-0.669437(5) 0.3070(3) 


SE 


-0.6693(1) 


0.307(1) 



" Based on n = {2, 4, 6} 
^ Based on n = {3, 5, 7} 
QMC (Quantum Monte Carlo) for square lattice^ 
SE (Series Expansion) for square lattice^ 

results. Thus in Table IIIII we show data for the case of the spin-1/2 HAF on the square 
lattice (corresponding to the case k = of the present Union Jack model). We present our 
CCM results in various LSUBn approximations (with 2 < < 7) based on the Union Jack 
geometry using the Neel model state. Results are given for the gs energy per spin E/N, 
and the magnetic order parameter M. We also display our extrapolated {n 00) results 
using the schemes of Eqs. ([6]) and (I7j) with the data sets n = {2,4,6} and n = {3,5,7}. 
The results are clearly seen to be robust and consistent, and for comparison purposes we 
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(a)n = {2,4,6} (b)n = {3,5,7} 

FIG. 4: (Color online) Ground-state energy per spin versus J2 for the Neel and canted 
phases of the spin-1/2 Union Jack Hamiltonian of Eq. ([1]) with Ji = 1. The CCM results 
using the canted model state are shown for various LSUBn approximations with (a) 
{n = {2, 4, 6} and (b) {n = {3, 5, 7}) with the canting angle (p = (pLsvBn that minimizes 
EhsvBn{(p)- We also show the n — 00 extrapolated result from using Eq. ([6]). 

also show the corresponding results using a QMC technique^ and from a linked-cluster series 
expansion (SE) method.- We note that for the square-lattice HAF no dynamic (or geometric) 
frustration exists and the Marshall- Peierls sign rule^ applies and may be used to circumvent 
the QMC "minus-sign problem". The QMC results^ are thus extremely accurate for this 
limiting (k = 0) case only, and represent the best available results in this case. Our own 
extrapolated results are in good agreement with these QMC benchmark results, as found 
previously (see, e.g., Ref. j38| and references cited therein) for CCM calculations performed 
specifically using the square lattice geometry, as well as for other CCM calculations for which 
the square-lattice HAF is a limit, such as for the spin-1/2 interpolating square-triangle J1-J2 
model,— for which the triangular lattice geometry was employed. It is gratifying to note in 
particular that although the individual LSUBn results for the spin-1/2 square-lattice HAF 
depend upon which geometry is used to define the configurations, the corresponding LSUBoo 
extrapolations are in excellent agreement with one another. 

In Fig. Hlwe show the CCM results for the gs energy per spin in various LSUBn approxi- 
mations based on the canted (and Neel) model states, with the canting angle 0LSUBn chosen 
to minimize the energy -E'LsuBn(0); as shown in Fig. [31 We also show separately the extrap- 
olated (LSUBoo) results obtained from Eq. ([6]) using the separate data sets n = {2,4,6} 
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FIG. 5: (Color online) Ground-state magnetic order parameter (i.e., the average on-site 

magnetization) versus J2 for the Neel and canted phases of the spin-1/2 Union Jack 
Hamiltonian of Eq. ([1]) with Ji = 1. The CCM results using the canted model state are 
shown for various LSUBn approximations with (a) [n = {2, 4, 6} and (b) n = {3, 5, 7}) 
with the canting angle (p = 4>hsvBn that minimizes -^LsuBn (</>)• We also show the n — )■ 00 

extrapolated result from using Eq. ([7]). 



and n = {3, 5, 7} as shown. As is expected from our previous discussion the energy curves 
themselves show very little evidence of the phase transition at k = ? with the energy and 
its first derivative seemingly continuous. 

Much clearer evidence for the transition between the Neel and canted phases is observed 
in our corresponding results for the gs magnetic order parameter M (the average on-site 
magnetization) shown in Fig. |5l For the raw LSUBra data we display the results for the Neel 
phase only for values of /t < k];^^^"' for clarity. However, the extrapolated (LSUBoo) results 
for the Neel phase are shown for all values of n in the range shown and where M > 0, using 
the extrapolation scheme of Eq. ([7]) and the LSUBn results based on the Neel model state. 
Once again we extrapolate the odd-n and even-n LSUBn results separately. It is interesting 
to note from Fig. |5]that the raw LSUBn data show the transition at ^LSUBn t^qj-q dearly for 
even values of n than for odd values of n. For the canted phase (for which 0LSUBn 7^ 0) we can 
clearly only show the extrapolated (LSUBoo) results using Eq. ([7]), for regions of k, for which 
we have data for all of the set n = {2, 4, 6} or n = {3, 5, 7}. We see from Table HT] that for the 
even-ra values we are limited (by the LSUB2 results) to values k, > k,^^^^"^ ^ 0.740, whereas 
for the odd-77, values we are limited (by the LSUB7 results) to values k, > k^^^^^ ~ 0.645. 
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The separate odd-n LSUBoo extrapolation curves for M for the Neel and canted phases are 
seen from Fig. [5b] to be extremely close at the value n = 0.645 and the curves appear to 
be about to meet at an angle which is either zero or very close to zero. A straightforward 
LSUBoo extrapolation using Eq. ([7]) for the three whole LSUBra curves (Neel plus canted) 
with n = {3, 5, 7} shows a value ~ 0.637 at which the extrapolated curve diverges (at 
zero or very small angle) from the corresponding LSUBoo estimate for the Neel state shown 
in Fig. |5bl For the corresponding even-ra data shown in Fig. [5al simple extrapolations of 
the LSUBoo curve to lower values of k < k,^^^^"^ ^ 0.740 using simple cubic or higher-order 
polynomial fits in n give a corresponding estimate of ~ 0.680 at which the Neel and 
canted phases meet. Both the even-n and odd-n LSUBn extrapolations yield a nonzero 
value for the average on-site magnetization of M ^ 0.195 ± 0.005 at the phase transition 
point ■ Thus the evidence from the behavior of the order parameter is that the transition 
at first-order one, in the sense that the order parameter does not go to zero at k^j, 

although it is certainly continuous at this point, and with every indication that its derivative 
as a function of k, is also continuous (or very nearly so) at k = k^. 

We also show in Fig. |6]the corresponding extrapolated (LSUBoo) results for the average 
on-site magnetization as a function of J2 (with Ji = 1), or hence equivalently as a function 
of K, for both the A sites (Ma) and the B sites (Mb) of the Union Jack lattice. We recall 
that, as shown in Fig. [Tal each of the A and B sites is connected to 4 NN sites on the square 
lattice by Ji-bonds, whereas each of the A sites is additionally connected to 4 NNN sites on 
the square lattice by J2-bonds. The extrapolations are shown in exactly the same regions, 
and for the same reasons, as those shown in Fig. [5l 

In Fig. [7] we also show the total gs magnetization per site, m = |(Mb — M^cos0), from 
using our CCM LSUBn results with = {2, 4, 6} and with the canting angle (p = ^LSUBn that 
minimizes -ELSUBn (0)- Clearly m = in the Neel phase where Ma = Mb and = 0. We also 
show in Fig. [7]the corresponding classical result m^^ = |s(l — cos 0ci) = i[l ~ (2^)"^] in the 
canted phase (with s = |). A comparison of the extrapolated (LSUBoo) CCM curve with 
its classical counterpart shows very clearly that the quantum fluctutions for this spin-half 
Union Jack model modify the classical behavior only relatively modestly, providing further 
evidence to what we also noted earlier in relation to Fig. |3l 

We also comment briefly on the large- J2 behaviour of our results for the canted phase. 
(We note that for computational purposes it is easier to re-scale the original Hamiltonian of 
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FIG. 6: (Color online) Extrapolated curves (LSUBoo) for the ground-state magnetic order 
parameters (i.e., the on-site magnetizations) Ma at sites A (joined by 8 bonds to other 
sites) and Mb at sites B (joined by 4 bonds to other sites) of the Union Jack lattice [and 

see Fig. [Ta] versus J2 for the Neel and canted phases of the spin-1/2 Union Jack 
Hamiltonian of Eq. ([T]) with Ji = 1. The CCM results using the canted model state are 
shown for various LSUBn approximations (n = {2, 4, 6} and n = {3, 5, 7}) with the canting 

angle = 0LsuBn that minimizes Ei^svBn{(f>) ■ 



Eq. ([T]) by putting J2 = 1 and considering small values of Ji.) The most interesting feature of 
the CCM results using the canted state as model state is that in all LSUBn approximations 
with n > 2 a termination point k^^^" is reached, beyond which no real solution can be 
found, very similar to the termination points shown in Fig. O For the even-n sequence the 
values are k^^^^ ^ 80 and k^^^^ ^ 80, whereas for the odd-ra sequence the values are 
^LSUBs ^ 250, kP"^^ ^ 85, and k^"^^ ^ 55. This is a first indication that the canted 
state becomes unstable at very large values of k against the formation of another (as yet 
unknown) state, as we discuss further in Sec. IIVBI below. 

Extrapolations of the gs energy using the data before the terminations points /tpuBn 
show that at large J2 values we have E/N — > — 0.3349J2 using the even-n LSUBn series and 
E/N — — O.3352J2 using the odd-n series. These numerical coefficients are precisely half of 
the values quoted in Table UTTl for the case J2 = 0. This is exactly as expected since both the 
K and the k — ?■ 00 limits of the Union Jack model are the square-lattice HAF, where 
in the latter case the square lattice contains only half the original sites, namely the A sites. 
Similarly, the extrapolated LSUBcx) values at larger values of k, before the termination point 
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FIG. 7: (Color online) The total ground-state magnetization per site, m = Mb — Ma cos 0, 
of the Union Jack lattice versus J2 of the spin-1/2 Union Jack Hamiltonian of Eq. ([1]) with 
Ji = 1. The CCM results using the canted model state are shown for various LSUBn 
approximations {n = {2,4,6}) with the canting angle = 0LSUBn that minimizes 
Ei,svBn{4>)- We also show the n — )■ 00 extrapolated result from using Eq. ((Tj) and compare 



for the on-site magnetization on the A sites are Ma — )■ 0.317 for the even-n LSUBn series 
and Ma — ?■ 0.306 for the odd-n series. Both values are again remarkably consistent with 
those shown in Table UTTl for the J2 = limit. The corresponding asymptotic values for the 
B-site magnetization are consistent with Mb — ?■ 0.5, as expected for large values of J2. 

B. Canted state versus the semi-striped state 

We turn finally to our CCM results based on the use of the semi-striped state shown in 
Fig. [lb] as the model state. Unlike in the case of the corresponding use of the canted state 
as model state, the results based on the semi-striped state do not terminate at a high value 
of J2 (with Ji = 1). We found no indication of such a termination value at any LSUBn 
level of approximation for 2 < n < 7 for values of J2 < 1000. All indications are thus that 
the semi-striped state is stable out to the J2 ^ 00 limit. Indeed for these LSUBn levels the 
CCM solutions based on the semi-striped state as model state exist for all values J2 > 1- 
For example, the LSUB6 solution based on the semi-striped state terminates at a lower 
end-point J2 ~ 0.41. 



it with the classical value m' 
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FIG. 8: (Color online) Comparison of the extrapolated (LSUBoo) curves for the 
ground-state energy per spin of the spin-1/2 Union Jack model, using the LSUBn data 
with n = {2, 4, 6} fitted to Eq. (|6]), for the CCM based on the canted and semi-striped 
states as model states. The LSUB4 and LSUB6 approximations (and hence also the 
extrapolated curve) for the canted model state terminate at J2 ~ 80. 

In Fig. Owe compare the extrapolated (LSUBoo) values of the gs energy per spin, based 
in each case on the LSUBtt, results with n = {2, 4, 6}, for our CCM results using the canted 
and semi-striped states. Although results for the canted model state become unavailable for 
J2 ^ 80 for the LSUB4 and LSUB6 approximations, the results based on the canted state lie 
lower in energy than those based on the semi-striped state for all values of J2 ^ 80 for which 
both sets of solutions exist. Although this is disappointing at first sight, the two sets of curves 
become extremely close for larger values of J2 as can be seen from Fig. |H1 Furthermore, we 
have also attempted a simple power-law extrapolation of the quantity E/{NJ2) for the gs 
energy of the canted state in powers of 1/^2? beyond the large- J2 LSUBn termination points 
(viz. at J2 ~ 80 for the LSUB6 approximation). Fits to sixth-, seventh-, and eighth-order 
polynomials give virtually identical results for values of J2 in the range 80 < J2 < 500 and 
these extrapolated curves do indicate that there is a second phase transition at ~ 125 ±5 
between the canted and semi-striped phases, such that for values k > the semi-striped 
phase becomes lower in energy. 

The gs energy in both the canted and semi-striped phases approaches the asymptotic 
value E/N ^ — 0.3349J2 for large values of J2 as J2 — )■ 00 (with Ji = 1). The corresponding 
asymptotic ( J2 — )■ 00) values for the average on-site magnetization of the semi-striped state 
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are Ma — )■ 0.317 for the A sites and Mb 0.5 for the B sites. 



V. DISCUSSION AND CONCLUSIONS 

In this paper we have used the CCM to study the influence of quantum fluctuations on 
the zero-temperature gs phase diagram of a frustrated spin-half Heisenberg antiferromagnet 
(HAF) defined on the 2D Union Jack lattice. We have studied the case where the NN Ji 
bonds are antiferromagnetic ( Ji > 0) and the competing NNN J2 = kJi bonds in the Union 
Jack array have a strength in the range < k < 00. On the underlying bipartite square 
lattice there are thus two types of sites, viz. the A sites that are connected to the 4 NN sites 
on the B-sublattice with Ji- bonds and to the 4 NNN on the A-sublattice with J2-bonds, and 
the B sites that are connected only to the 4 NN sites on the A-sublattice with Ji-bonds. The 
K, = limit of the model thus corresponds to the spin-half HAF on the original square lattice 
(of A and B sites), while the k — 00 limit corresponds to the spin-half HAF on the square 
lattice comprised of only A sites. We have seen that at the classical level (corresponding to 
the case where the spin quantum number s — )■ 00) this Union Jack model has only two stable 
gs phases, one with Neel order for n < = 0.5 and another with canted ferrimagnetic order 
for K, > nf. We have therefore first used these two classical states as CCM model states to 
investigate the effects of quantum fluctuations on them. 

For the spin-half model we find that the phase transition between the Neel antiferromag- 
netic phase and the canted ferrimagnetic phase occurs at the higher value = 0.66 ±0.02. 
The evidence from our calculations is that the transition at subtle one. From the 

energies of the two phases it appears that the transition is either second-order, as in the 
classical case, or possibly, weakly first-order. However, on neither side of the transition at 
Kci does the order parameter M (i.e., the average on-site magnetization) go to zero. Instead 

clS Ki y 

from either side, M — )■ 0.195 ± 0.005, which is more indicative of a first-order 
transition. Furthermore, the slope dM / dn of the average on-site magnetization as a function 
of K also seems to be either continuous or to have only a very weak discontinuity at k = • 
Before continuing with the possibility of a further phase we compare our results with 
those from previous calculations of the same model using spin-wave theory (SWT)^'2I and 
the linked-cluster series expansion (SE) method.— Collins et al.—^ used linear (or leading- 
order) spin-wave theory (LSWT) to show that on the basis of a comparison of the gs energies 
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of the two phases, the phase transition between the Neel and canted phases is of first-order 
type and occurs at ~ 0.84. In LSWT the Neel staggered magnetization per site M 
remains substantial at this estimate of k^^ ^ 0.84. However, it is well known that LSWT 
results become unreliable near the transition region, and they surmised that the Neel order 
parameter M might vanish at or before this point, yielding a possible scenario where a 
second-order Neel transition might occur at a value ^ 0.84, followed by a possible 
intermediate spin-liquid phase (as in the pure J1-J2 model, as discussed in Sec. [T]), and 
then a first-order transition to the canted phase at a somewhat larger value of k. Our own 
results provide no evidence at all for such an intermediate spin-liquid phase between the 
Neel antiferromagnetic and canted ferrimagnetic phases. 

Although LSWT is known to give a reasonable description of the spin-1/2 Heisenberg 
antiferromagnet on the square lattice (k = 0), it is surely unable to model the frustrated or 
intermediate regime accurately. Similar shortcomings of spin-wave theory (SWT) have been 
noted by Igarashi^ in the context of the related spin-1/2 J1-J2 model on the square lattice, 
discussed briefly in Sec. [B He showed that whereas its lowest-order version (LSWT) works 
well when J2 = 0, it consistently overestimates the quantum fluctuations as the frustration 
J2/ Ji increases. In particular he showed, by going to higher orders in SWT in powers of 1/s 
where s is the spin quantum number and LSWT is the leading order, that the expansion 
converges reasonably well for J2/ Ji < 0.35, but for larger values of J2/ Ji, including the point 
J2/ Ji = 0.5 of maximum classical frustration, the series loses stability. He also showed that 
the higher-order corrections to LSWT for J2/J1 ^ 0.4 make the Neel-ordered phase more 
stable than predicted by LSWT. He concluded that any predictions from SWT for the spin- 
1/2 J1-J2 model on the square lattice are likely to be unreliable for values J2/ Ji ^ 0.4. It is 
likely that a similar analysis of the SWT results for the spin-1/2 Union Jack model studied 
here would reveal similar shortcomings of LSWT as the frustration parameter k = J2/J1 is 
increased. 

In a later paper— by Zheng et al, SE techniques were applied to our spin-half Union Jack 
model and were compared with those from both LSWT for both the Neel and canted phases 
and modified second-order SWT for the Neel phase. Using the SE method for the Neel 
phase gave what these authors termed very clear evidence of a second-order phase transition 
at a critical coupling = 0.65 ± 0.01 at which the Neel staggered magnetization per site 
vanished. For higher couplings the system was seen to lie in the canted phase, with no sign 
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of any intermediate spin-liquid phase between these two magnetically ordered states. Use 
of the SE method in the canted phase produced a gs energy which continues smoothly from 
the Neel into the canted phase. Zheng el ai— found, furthermore, that in the canted phase 
the staggered magnetizations per site in both the vertical and horizontal directions shown in 
Fig. [Talalso appear to drop smoothly toward zero around the same value = 0.65 ± 0.01, 
albeit with very large error bars. 

The above SE estimate for is clearly in excellent agreement with our own. However, 
whereas the evidence from the order parameter M from the SE technique clearly favors a 
second-order transition at at which M — )■ from both sides, our own CCM calculations 
clearly favor a first-order transition at which M — )■ 0.195 ± 0.005. We note, however, that 
the errors on the SE estimates for M become increasingly large as the phase transition at 
Kci is approached from either side. We believe that this could easily account for the seeming 
discrepancy between our respective predictions for the order of the phase transition at Kci- 
We note too that Zheng el alr^ were themselves puzzled by the discrepancy between the 
prediction of SWT that the Neel magnetization per site M does not vanish at and that 
of the SE technique that M vanish does vanish there. While they recognized (as do we, as 
we discussed above) that SWT cannot be taken as an infallible quide, they found the huge 
difference with the prediction from the SE technique perturbing. Those authors ended by 
stating that, in their opinion, the nature of the transition from the Neel to the canted phase 
in the spin-half Union Jack model deserved further exploration. We believe that our own 
work reported here has considerably illuminated the transition at Kci- 

Neither SWT nor SE techniques have been applied to the possible semi-striped state of 
Fig. [lb] for the spin-half Union Jack model and so we have no results against which to com- 
pare our own. We were led to consider such a state as a possible gs phase of the model at 
large values of k as discussed in Sec. [Tll Thus, to recapitulate, the k — )■ oo limit of the canted 
phase of the Union Jack model (for either the quantum s = 1/2 model considered here or the 
classical s — )■ oo case) gives a state in which the spins on the antiferromagnetically-ordered 
A-sublattice are orientated at 90° to those on the ferromagnetically-ordered B-sublattice. 
The actual k — oo limit should, in either case, be decoupled antiferromagnetic (A) and ferro- 
magnetic (B) sublattices, with complete degeneracy at the classical level for all angles of rel- 
ative ordering directions between the two sublattices. We argued that quantum fluctuations 
could, in principle, lift this degeneracy by the well-known order by disorder phenomenon.— 
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Since quantum fluctuations are also well known from many spin-lattice problems to favor 
collinearity, there is a strong a priori possibility that the semi-striped state of Fig. [lb] might 
be energetically favored at large values of k over the non-coUinear state which is the k ^ oo 
limit of the canted state in which — )■ 90°. 

Accordingly we repeated our CCM calculations using the semi-striped state as model 
state. We found some evidence that at very large values of k there might indeed be a second 
phase transition at ~ 125±5, based on the relative energies of the canted and semi-striped 
states. Such a prediction is based, however, on an extrapolation of the data on the canted 
state into regimes where the CCM equations have no solution for LSUBn approximations 
with n > 3, and hence cannot be regarded as being as reliable as our prediction for Kci- 
If the phase transition at k,c2 does exist it would be of first-order type according to our 
results. It would be of considerable interest to explore the possible transition at k,c2 between 
the canted and semi-striped phases by other techniques, possibly including SWT and SE 
methods. 

As has been noted elsewhere,— high-order CCM results of the sort presented here have 
been seen to provide accurate and reliable results for a wide range of such highly frustrated 
spin-lattice models. Many previous applications of the CCM to unfrustrated spin models 
have given excellent quantitative agreement with other numerical methods (including exact 
diagonalization (ED) of small lattices, quantum Monte Carlo (QMC), and series expansion 
(SE) techniques). A typical example is the spin-half HAF on the square lattice, which is the 
K = limit of the present model (and see Table [In|) . It is interesting to compare for this k = 
case, where comparison can be made with QMC results, the present CCM extrapolations 
of the LSUBn data for the infinite lattice to the n — )■ oo limit and the corresponding QMC 
or ED extrapolations for the results obtained for finite lattices containing spins that 
have to be carried out to give the N ^ oo limit. Thus, for the spin-1/2 HAF on the square 
lattice the "distance" between the CCM results for the ground-state energy per spin^ at the 
LSUB6 (LSUB7) level and the extrapolated LSUBoo value is approximately the same as the 
distance of the corresponding QMC result^ for a lattice of size = 12 x 12 (A^ = 16 x 16) 
from its A^ — oo limit. The corresponding comparison for the magnetic order parameter 
M is even more striking. Thus even the CCM LSUB6 result for M is closer to the LSUBcxo 
hmit than any of the QMC results for M for lattices of A^ spins are to their A^ — )■ oo limit 
for all lattices up to size A^ = 16 x 16, the largest for which calculations were undertaken.— 
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Such comparisons show, for example, that even though the "distance" between our LSUBn 
data points for M and the extrapolated [n — )■ oo) LSUBoo result shown in Fig. [5] may, at 
first sight, appear to be large, they are completely comparable to or smaller than those in 
alternative methods (where those other methods can be applied). Furthermore, where such 
alternative methods can be applied, as for the spin-1/2 HAF on the square lattice, the CCM 
results are in complete agreement with them. 

By contrast, for frustrated spin- lattice models in two dimensions both the QMC and ED 
techniques face formidable difficulties. These arise in the former case due to the "minus-sign 
problem" present for frustrated systems when the nodal structure of the gs wave function is 
unknown, and in the latter case due to the practical restriction to relatively small lattices 
imposed by computational limits. The latter problem is exacerbated for incommensurate 
phases, and is compounded due to the large (and essentially uncontrolled) variation of the 
results with respect to the different possible shapes of clusters of a given size. 

For highly frustrated spin-lattice models like the present Union Jack model, a powerful 
numerical method, complementary to the CCM, is the linked-cluster series expansion (SE) 
technique.—"—'^"— The SE technique has also been applied to the present model.— Our own 
results have shed considerable light on the nature of the phase transition at hi^i observed by 
SE techniques and the discrepancies between the results from SE and SWT methods. 

We end by remarking that it would also be of interest to repeat the present study for the 
case of the s > 1/2 Union Jack model. The calculations for this case are more demanding 
due to an increase at a given LSUBn level of approximation in the number of fundamental 
configurations retained in the CCM correlation operators. Nevertheless, we hope to be able 
to report results for this system in the future. 
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